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CLOSE-PACKED, UNIFORMLY ADJACENT, MULTERESOLUTIONAL, 
OVERLAPPING SPATIAL DATA ORDERING 

Description 
FIELD OF THE INVENTION 

5 The invention refers generally to methods and digital information systems that represent 

visual images of planar or near planar space. More particularly, the methods and systems 
pertain to images that are referenced to earth. 
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Description of Prior Art 

Storing Images as Discrete Grids 

A common method of representing spatial information in a digital device is the division 
30 of space into a tessellation of discrete cells. Examples include raster images and pixel based 
visual devices. Raster images are visual references to data in the form of "dots"-typically 
arranged in an array of rows and columns, along with an attribute of that dot-usually a number 
representing a sensor's spectral value. 
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Vector images are interpretations of spatial features using geometric shapes described 
mathematically with coordinates. Identification of point locations, the basic unit of vector 
mathematics, is often based on the rectangular Cartesian grid. 

Regular rectangular grids are commonly used due to their immediate transformation into 
5 simplified matrix or array operations, their compatibility with well-understood planar Euclidean 
geometry, and the prevalence of cardinal directions in memetic spatial thought. 

Systems that utilize Spatial Data Representation 

A Geographic Information System [GIS] is a computer application that provides for the 
creation, storage, visualization and analysis of gee-referenced spatial information. A GIS can be 

10 raster or vector based. In a GIS, spatial references are made to geographic locations on the earth 
surface. This reference location is further enhanced with descriptive information regarding the 
feature or area it represents. Georeferenced spatial information encompasses elements familiar 
to those using paper maps-such as roads and contours and place names, but also includes any 
data with a physical location description-like street addresses, cadastral boundaries and 

15 geographic coordinates and can include space and airborne imagery. Almost all data can 
contain geospatial elements. 

Some Geographic Information Systems incorporate 3-dimensions wim a Digital 
Elevation Model, a specialized form of geospatial vector data that describes the 3-dimensions of 
a terrain and heights on features. A Discrete Global Grid System is a type of GIS that provides 

20 a tessellation of uniform cells placed over the earth as a reference to location and area of 
influence. 

Hexagonal Grids 

Hexagonal cells are documented as the preferred shape for many gridding applications. 
A hexagonal grid is the most efficient method of dividing planar space into regions of equal area 

25 with the least total perimeter. For a given distance between the center of a polygon and its 
farthest perimeter points, the hexagon has the largest area of the three possible uniform cell 
grids, rectangle, triangle and hexagon. The hexagonal grid provides a symmetric cell structure 
such that all adjacent cells are equidistant to 6 edge neighbours. Hexagons provide the peak 
efficiency for spatial data storage. Hexagon lattices have been shown to provide excellent 

30 modeling characteristics. 
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Indexing systems have been proposed for transposing spatial imagery on hexagon grids, 
notably Generalized Balanced Ternary of Gibson and Lucas. Hexagons can be used to define 
geographic location. Sahr, et. al. proposed a Discrete Global Grid System called the 
Icosahedron Snyder Equal Area Aperture 3 Hexagon Grid as an alternative to rectangular 
5 graticule of latitudes and longitudes based on a projection from the icosahedron to the spheroid 
proposed by John Snyder. 

Hexagonal grids have been overlooked in many applications due to their 
nonconformance, by definition, to rectangular gridding and the utility of accompanying 
rectangular mathematics. Hexagonal Grids have presented an inability to aggregate completely 
10 into larger, or disaggregate completely into smaller, self-similar cells in hierarchal B-Tree 
divisions such as the rectangular Quad Tree or sub-forms, such as the Quaternary Triangular 
Mesh, proposed by Geoffrey Dutton. 

Addressing with Hierarchal Indices 

On-Line Analytical Processing is a term used to define systems that promote fast 
15 analysis of shared multidimensional information. Prior art, in the development of relational and 
spatial database, place heavy reliance on balanced hierarchal indexing for efficient data search 
and retrieval. In regular graticulation, the division of the earth surface in latitudes and 
longitudinal angles, hierarchies are not apparent nor are the resulting cells of equal area or 
shape. 

20 Spatial Transformations 

Spatial analysis, data fusion and modeling represent sophistication in information 
systems classified as transformations. In a GIS, spatial analysis is generally based on ideas of 
proximity, networks and neighbours along with standard data analysis techniques such as 
correlative analysis. Analytical transforms tend to cut the data in several ways to facilitate 

25 interpretation. 

Data fusion has many definitions. Generally, it can be held to mean the creation or 
enhancement of georeferenced attributes based on the combination of two or more attributes. 
An example would be the use of high-resolution panchromatic satellite imagery to enhance the 
values of low-resolution hyper-spectral satellite imagery. 
30 Spatial modeling provides a means of predicting outcomes based on known behaviors 

and the presence of specified factors. Modeling uses known behaviors to extrapolate or 
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simulate possible outcomes. Finite difference methods can utilize standardized meshes to route 
partial differential equations to predict weather patterns. Cellular automata require simple rules 
and cell boundaries to predict future behaviors and patterns as demonstrated by John Conway in 
his cellular automata discovery: "the game of life". A practical use the hexagon grid as a 
5 framework to complete environmental systems modeling is introduced by U.S. Panu in a semi- 
deterministic hexagonal based runoff model. 

Stephen Wolfram has shown the successful use of hexagonal meshes for using cellular 
automata to simulate systems described by partial differential equations such as those that 
describe the flow of fluid, diffusion or heat transfer. 

10 On Demand On-line Data Transmission 

Information systems are ideally utilized when the data is available to the maximum 
number of users. Earth observation systems proliferate global data sets. On-line on-demand 
data systems that promote file sharing can be used in a spatial data environment. Means of 
storage, retrieval, searching and finding the data is necessary to maximize data utility. Peer-to- 

15 peer file sharing and data discovery applications are well known in prior art. Such systems 
require a client to search repositories of data describing the address of server and content of data 
files. When the search criteria are met, the address of the server where the data is stored is 
returned. 

Since these data sets can be large, development of data compression, peer-to-peer, and 
20 streaming data file sharing and grid computing applications that provide opportunities for . 
efficient file transfer and shared processing are beneficial. Progressive transmission and 
continuous transmission are useful ways of selectively transmitting data. 

Current Situation and Present Need 

The planet earth is a dynamic sphere containing interrelationships of environmental 
25 systems: an arrangement of geophysical, meteorological and biological activity. Within these 
natural systems, we create divisions of human interest: artificial boundaries of geopolitical, 
functional or demographic areas, and natural boundaries such as watersheds or resource 
delineations. Also, human assets or infrastructure like roads, buildings, sewers and landfill 
sites, serve as detailed areas of interest within human biological habitat. 

30 Many areas of spatial understanding such as epidemiology, resource management and 

biodiversity would improve with the availability of detailed geospatial data at multiple 
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resolutions coupled with the behavioral understandings of stochastic methods to provide semi- 
deterministic modeling. 

As space, airborne and ground-based sensors increase in number and sophistication; 
there is a rapid, and potentially unlimited, growth of global information characterizing these 
5 environmental systems. The GIS-centric approach to geo-spatial data and applications as a 
paradigm of static well defined spatial data is being eclipsed by these realities and the current 
operational needs of many spatial data users. A re-conceptualization of the way geo-spatial 
information is processed toward a general case where all geo-spatial is dynamic and real-time 
will necessitate a corresponding effective and practical way of handling these diverse and 
10 increasingly large, spatial data sets. 

Typically Geographic Information Systems reference data to flat maps by projecting 
geographic coordinates to a regular rectangular grid. The distortion that results from 
referencing the geodesic sphere onto a projected regular grid map, like the Mercator Projection, 
is a characteristic that school children are familiar with. A similar condition produces 
15 significant difficulties in global geospatial data storage. Traditional geographic coordinate 
systems, latitude and longitude, do not form uniform cells and therefore are not idealized for use 
where uniformity of cells is assumed or is beneficial. 

Merging geospatial information of disparate sources, accuracies, precisions, and formats 
is complex. Unlike textual or fielded data, spatial temporal data requires a common spatial 
20 reference and significant background or meta-data to facilitate cut and paste, overlay and mosaic 
operations. Solving the problem of integrating data requires an idealized common reference 
model. As raster, vector and digital elevation models rely on a regular grid, a standard reference 
grid would serve as a framework for the compilation and enhancement of all geospatial 
information. 

25 To access global data, an idealized model of the earth would best approximate it, as it is, 

an irregular spheroid. This idealized model would allow characteristics of the earth at any point 
in time and at any location to be observed. An infinite resolution of uniform discrete cells 
provides a mechanism to build such a data model. A digital version of the globe as a 
visualization application could show the world in a space view and, as the user rotates the globe 

30 and zooms in to a specific point of interest, the spatial data that it represents would confer the 
curvature of the surface onto the computer monitor in increasing detail. 
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SUMMARY OF THE INVENTION 

It is therefore an object of the invention to provide for an efficient method of mapping 
two-dimensional spatially organized data to the one-dimensional linear space of a computer 
storage medium. This object is solved by a method according to claim 1 and a discrete global 
grid system according to claim 7, claims 2 to 6 refer to specifically preferred realizations of the 
inventive method, claims 8 to 21 refer to specific advantageous realizations of the discrete 
global grid system according to claim 7, whereas by the features of claim 14 especially a spatial 
data retrieval system can be realized, by the features of claim 15 a digital globe visualization 
system can be realized, by the features of claim 16 an online spatial data discovery and file 
sharing tool can be realized, by features of claim 17 a spatial data browser can be realized, by 
the features of claim 18 a spatial data analyzer can be realized and by the features of claim 19 a 
geospatial model building system can be realized. 

The invention realizes an efficient method of mapping two-dimensional spatially 
organized data to the one-dimensional linear space of a computer storage medium using an 
15 assignment of spatial attributes to a cellular tessellation. A linear indexing provides a 
hierarchical data structure for closely packed, uniformly adjacent, multiresolutional, overlapping 
cells. The method promotes the hexagonal grid as a practical alternative to rectangular or 
triangular grid use in a digital computer environment. 

According to one aspect of the invention is a method for storing two-dimensional 
20 spatially organized data in one-dimensional space on a computer storage medium by mapping 
the attributes of continuous state planar space to a multi-resolutional tessellation of closely 
packed uniform cells, each cell being uniquely identified with a sequential number whereas the 
number includes the identification of a parent cell, the parent cell encompassing a cluster of 
child cells in a spatial hierarchy of specific order thereby identification of neighbour cells and 
25 child cells comprising the requirements: 

Req. 1. Spatial attributes are assigned to a parent cell, whose centroid represents its location and 
the voronoi region created by the boundary with adjacent parent centroids forming the 
closed area for which the properties of the cell are represented; 
Req. 2. A parent cell for which the centroid location is not a centroid location for any lower 
30 resolution cells defines the location of a single new child cell of the next highest 

resolution; alternatively, 
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Req. 3. A parent cell for which its centroid location is also a centroid location for any lower 
resolution cells defines the location of a single new child cell of Ihe next highest 
resolution and multiple new child cells of the next highest resolution, one located at each 
of the vertices of the parent's boundary edge. 
5 Whereby during initial conditions, a parent cell will be assigned a general hexagon shape 

or the shape of the plane for which it represents, with a starting centroid location that can be 
considered the planar origin. 

Embodiments of the invention include an on-line Discrete Global Grid System featuring 
hexagonally gridded data for efficient computational storage, retrieval, integration, transmission, 
10 visual display, analysis, fusion, and modeling. 

BRIEF DESCRIPTION OF THE DRAWINGS 
FIG. 1 - TYPICAL COMPUTER HARDWARE FOR PYXIS OPERATION - illustrates an 
apparatus and operating environment in conjunction with which embodiments of the invention 
may be practiced. 

15 FIG. 2 - PRIOR ART FORM OF SPATIAL DATA ON A RECTANGULAR PLANE - 
illustrates prior art methods of locating spatial information on a rectangular grid. 
FIG. 3 - SPATIAL DATA ON A PYXIS HEXAGON PLANE - Illustrates an alternative planar 
location scheme called the PYXIS plane whereas the points are closely packed. 
FIG. 4a, 4b, 4c, 4d, 4e and 4f - EXEMPLORY ORDERING OF THE CHILD CLUSTERED 

20 CELLS 

FIG. 5a, 5b, 5c, 5d, 5e and 5f - SIX SUCCESSIVE RESOLUTIONS OF PYXIS INDEXING 
USING SIMPLE ORDERING 

FIG. 6a and 6b - PYXIS HIERARCHY - illustrates hierarchal views of the PYXIS indexing. 
FIG. 7 - TRANSLATIONAL MOVEMENT AND ADDITION ON PYXIS PLANE - illustrates 
25 an exemplary procedure with addition table for identifying nearest neighbour cells and other 
such operations on the PYXIS plane. 

FIG. 8 - BOOLEAN OPERATIVES ON PYXIS PLANE - illustrates an exemplary procedure 
for identifying a PYXIS point as spatial contained, overlapping or excluded, from a specific cell 
area. Figures included are FIG. 8a - PYXIS Plane FIG. 8b- Successive Cell Resolutions, FIG. 
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8c. - Resolution Dependent Conditions within Overlapping Area on Vertex - Condition 2.1, 
FIG. 8d. - Resolution Dependent Conditions within Overlapping Area on Edge - Condition 2.2. 

FIG. 9a - PYXIS PLANE ON TRIANGULAR FACE, 9b - PYXIS PLANE ON 
ICOSAHEDRON, 9c - PENTAGONAL LAYOUT OF PYXIS PLANE, 9d - PYXIS PLANE 
5 ON UNFOLDED ICOSAHEDRON and 9e - PYXIS PLANE ON AUTHALIC SPHERE - 
illustrate an exemplary method for establishing and referencing multiresolutional closely packed 
cells as a digital reference and representation of the earth spheroid by projecting the PYXIS 
plane grid from a geodesic solid, the icosahedron. 

FIG. 10 - DISCRETE GLOBAL GRID SYSTEM - is an illustration of the components of a 
10 Geodesic Discrete Global Grid System. 

FIG. 11- EXAMPLE OF RECTANGULAR TO PYXIS CONVERSION - is an illustration of 
procedure for determining Rectangular Point containment in PYXIS Cell. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
Apparatus and Operating Environment 

15 FIG. 1 is a diagram of the apparatus and operating environment in conjunction with 

which embodiments of the invention may be practiced. The description of FIG. 1 is intended to 
provide a brief; general description of suitable computer apparatus in conjunction with which 
the invention may be implemented. The computer, in conjunction with which embodiments of 
the invention may be practiced, may be a conventional computer, a distributed computer, an 

20 embedded computer or any other type of computer; the invention is not so limited. Such a 
computer typically includes one or more processing units as its processor, and a computer- 
readable medium such as a memory. The computer may also include a communications device 
such as a network adapter or a modem, so that it is able to communicatively couple other 
computers. 

25 Although not required, the invention is described in the general context of computer- 

executable instructions being executed by a computer. Generally, program modules include 
routines, programs, objects, components, data structures, etc. that performs particular tasks or 
implement particular abstract data types. 

Moreover, those skilled in the art will appreciate that the invention may be practiced 
30 with other computer system configurations, including hand-held devices, multiprocessor 
systems, microprocessor-based or programmable consumer electronics, network PCS, 
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minicomputers, mainframe computers, and the like. The invention may also be practiced in 
distributed computing environments where tasks are performed by remote processing devices 
that are linked through a communications network. In a distributed computing environment, 
program modules may be located in bom local and remote memory storage devices. 
5 The exemplary hardware and operating environment of FIG. 1 for implementing the 

invention includes a general purpose computing device in the form of a computer (1), including 
a processing unit (2), a system memory (3), and a system bus (4) that operatively couples 
various system components, including the system memory (3), to the processing unit (2). There 
may be only one or mere may be more than one processing unit (2). 
10 The system bus (4) may be any of several types of bus structures including a memory 

bus or memory controller, a peripheral bus, and a local bus using any of a variety of bus 
architectures. The system memory (3) includes read only memory (5) and random access 
memory (6). A basic input/output system (BIOS) (7), containing the basic routines mat help to 
transfer information between elements within the computer (1), such as during start-up, is stored 
15 in random access memory (6). The computer (1) further includes a hard disk drive (8), a 
magnetic disk drive (9) and an optical disk drive (10), as memory storage devices. 

The hard disk drive (8), magnetic disk drive (9), and optical disk drive (10) are 
connected to the system bus (4) by a hard disk drive interface (11), a magnetic disk drive 
interface (12), and an optical disk drive interface (13), respectively. The drives and their 
20 associated computer-readable media provide nonvolatile storage of computer-readable 
instructions, data structures, program modules and other data for the computer (I). It should be 
appreciated by those skilled in the art that any type of computer-readable media which can store 
data that is accessible by a computer, such as magnetic cassettes, flash memory cards, digital 
video disks, Bernoulli cartridges, random access memories (RAMs), read only memories 
25 (ROMs), and the like, may be used in the exemplary operating environment. 

A user may enter commands and information into the personal computer (1) through 
input devices such as a keyboard (14) and mouse (15). Other input devices (not shown) may 
include a microphone, joystick, game pad, satellite dish, scanner, or the like. These and other 
input devices are often connected to the processing unit (2) through a serial port interface (16) 
30 that is coupled to the system bus, but may be connected by other interfaces, such as a parallel 
port, game port, or a universal serial bus (USB). A monitor (17) or other type of display device 
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is also connected to the system bus (4) via an interface, such as a video adapter (18). In addition 
to the monitor, computers typically include other peripheral output devices (not shown), such as 
speakers and printers. 

The computer (1) may operate in a networked environment using logical connections to 
5 one or more remote computers, such as remote computer (19). A communication device coupled 
to or a part of the computer (1) achieves these logical connections; the invention is not limited to 
a particular type of communications device. The remote computer (19) may be another 
computer, a server, a router, a network PC, a client, a peer device or other common network 
node, and typically includes many or all of the elements described above relative to the 
10 computer 1. The logical connections depicted in FIG. 1 include a local-area network (LAN) 
(20) and a wide-area network (WAN) (21). Such networking environments are commonplace in 
office networks, enterprise-wide computer networks, intranets and the Internet, which are all 
types of networks. 

When used in a LAN-networking environment, the computer (1) is connected to the 
15 local network (20) through a network interface or adapter (22), which is one type of 
communications device. When used in a WAN-networking environment, the computer 1 
typically includes a modem (23), a type of communications device, or any other type of 
communications device for establishing communications over the wide area network (21), such 
as the Internet. The modem (23) is connected to the system bus (4) via the serial port interface 
20 (16). In a networked environment, program modules depicted relative to the personal computer 
(1), or portions thereof, may be stored in the remote memory storage device. It is appreciated 
that the network connections shown are exemplary and other means of and communications 
devices for establishing a communications link between the computers may be used. 

A Method to Structure Spatial Data at titeform of a Linear Index 

25 Division of Space 

FIG. 2 illustrates prior art methodology for locating spatial information on a rectangular 
grid. A single spatial coordinate, Point (24), from a defined rectangular Plane (25), with 
Attributes (26), such as on, off, red, blue, hot, cold, 20, 30, 40, etc, shown. Point (24) and 
Attributes (26) have been created, collected or recorded by any of a number of analog and/or 

30 digital methods. Point (24) can be representative of any single dot or a singular dot within a 
large group of a Raster Image (27) located with a row and column. Alternatively, Point (24) 
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could be a point of a mathematically defined Vector Line (28) located with a Coordinate (6,7) 
(29) and an index Link (30) to Attributes (26). Exemplary linear and area constructs of the 
point include a Network (31) and Polygon (32). 

Point (24), by its nature, also represents a centroid and a cell area that encompasses a 
5 distance between it and the next adjacent point or, on the basis of known practice within 
geometry and other mathematics, Point (24) is bound by its accuracy, error and precision based 
on factors of its creation and record. It is exemplary that within this cell area, Attributes (26) 
can be considered homogenous. For the purpose of this description, the area, whether implied 
explicitly or implicitly, will be referred to as a point's area of influence. 

10 FIG. 3 illustrates an alternative planar location scheme whereas the points are closely 

packed. Referring to FIG. 3, the single spatial Coordinate (29) of Point (24), is now defined by 
a hexagonal PYXIS Plane (33), as Cell (34) with the same Attributes (26), on, off, red, blue, hot, 
cold, 20, 30, 40, etc. 

Successive aperture three subdivisions of a general hexagon cell define locations on the 
15 PYXIS Plane (33). In an aperture three hexagon subdivision, at any given resolution "n", the 
area of a general hexagon is 3 times the area of a hexagon one resolution higher, "n+1". 

The Plane (33) is infinite in two-dimensional space. As also illustrated on FIG. 2, Cell 
(34) can also be representative of any single dot or a singular dot within a large group of a 
Raster Image (35) or Cell (34) could be a point of a Vector Line (28), either located with a 
20 PYXIS Index (7277472) (36), which is explained in more detail below, and used as the index to 
Attribute (26). The area bound by Cell (34) becomes its area of influence. 

The division of space on the PYXIS Plane (33), herein coined as the PYXIS innovation, 
adheres to the following requirements: 

Req. 1. Spatial attributes are assigned to a parent cell, whose centroid represents its location and 
25 the voronoi region created by the boundary with adjacent parent centroids forming the 

closed area for which the properties of the cell are represented; 

Req. 2. A parent cell for which the centroid location is not a centroid location for any lower 
resolution cells defines the location of a single new child cell of the next highest resolution; 
alternatively, 
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Req. 3. A parent cell for which its centroid location is also a centroid location for any lower 
resolution cells defines the location of a single new child cell of the next highest resolution 
and multiple new child cells of the next highest resolution, one located at each of the 
vertices of the parent's boundary edge. 
5 Whereby during initial conditions, a parent cell will be assigned a general hexagon shape 

or the shape of the plane for which it represents, with a starting centroid location that can be 
considered the planar origin. A note of interest is that the resulting outline is a figure with 
properties of the Koch Snowflake curve-amongst others, a infinite perimeter with finite area. 

Addressing and Indexing 

10 The system or procedure for addressing the cells is based on assigning a number 

consisting of the cell index of its related parent plus the number corresponding to its order in the 
cluster of child cells. Ordering of the cluster of child cells can follow simple sequential ordering 
or other exemplary methods shown: Generalized Balanced Ternary FIG. 4a, Gray Coding FIG. 
4b, simple ordering FIG. 4c, various z curves FIG. 4d, 4e, and 4f and other forms known to 

15 those versed in the art such as modifications of Morton, Hilbert, Koch, Peano curves and/or 
ordering that are not shown. 

Child cells at a given resolution (n+1) are uniquely related to a next lower parent cell at 
resolution (n) if they share the same centroid location or, if the centroid of the child cell is 
located on the vertex of the cell boundary of a next lower parent cell at resolution (n) whose 

20 centroid location is also a centroid location of its parent cell(s) at any lower resolution (<n-l). 

FIG. 5a to 5f illustrate six successive resolutions of the PYXIS indexing using a simple 
ordering. 

FIG. 5a. A generalized hexagon is set at the origin as an initial condition with a given 
Resolution (37) and assigned a place value of 7. 

25 FIG. 5b. The second Resolution (38) cell locations exist at all centroids and vertices of 
Resolution (37) cells, following rule 3, and are indexed with the 7 of the parent 
and the ordering 1, 2, 3, 4, 5, 6 (of the vertices) and following 7 (at the centroid) 
resulting in new indices 71, 72, 73, 74, 75, 76 and 77. 

FIG. 5c. The third Resolution (39) cell locations follow the same pattern; new cells are 
30 given the value of the next lower Resolution (38) parent and a centroid value (7) 
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if they share the same centroid location, rule 2 (71, 72, 73, 74, 75, and 76 and 77 
spawn centroid cells 717, 727, 737, 747, 757, and 767), and if the centroid of the 
new cell is located on the vertex of the cell boundary of a Resolution (38) cell 
whose centroid location is also a centroid location for a Resolution (37) cell, rule 
5 3, the new cells are given the parent value and an order value (77 spawn vertex 

cells 771, 772, 773, 774, 775, 776 and 777). 

FIG. 5d. The forth Resolution (40) cell division continues similarly, all new cells are 
given the value of the next lower Resolution (39) index if they share the same 
centroid location, rule 2 (771, 772, 773, 774, 775, and 776 spawn centroid cells 

10 7717, 7727, 7737, 7747, 7757, and 7767), and if the centroid of a new cell is 

located on the vertex of the cell boundary of a Resolution (39) cell whose 
centroid location is also a centroid location for a Resolution (38) cell, rule 3, the 
new cells are given the parent value and an order value (717, 727, 737, 747, 757, 
767 and 777 spawn vertex cells 7171, 7172, 7173, 7174, 7175, 7176, 7177, 7271, 

15 7272, 7273, 7274, 7275, 7276, 7277, 7371, 7372, 7373, 7374, 7375, 7376, 7377, 

7471, 7472, 7473, 7474, 7475, 7476, 7477, 7571, 7572, 7573, 7574, 7575, 7576, 
7577, 7671, 7672, 7673, 7674, 7675, 7676, 7677, 7771, 7772, 7773, 7774, 7775, 
7776, and 7777). 

FIG. 5e. This division is continued with successively finer Resolution (41) and 
20 FIG. 5f. Resolution (42). 

It should be noted that to facilitate available space, indices for new cells located on 
centroids used in lower resolutions are not shown on FIG. 5d, 5e, and 5f. 

These unique indices provide for the mapping of two-dimensional spatially organized 
data in the form of multiple resolutions of closely packed uniformly adjacent cells to a linear 
25 index compatible for use in the one-dimensional linear space of a computer memory and/or 
storage medium, as a conventional file or relational, multidimensional or other similar database 
system familiar to those versatile in the art. Also, the indexing is so structured as to provide a 
common addressing scheme for traditional forms of spatial data with an explicit indication of 
the spatial data location and its area of influence. 

30 As spatial data structures require some flexibility, exceptions to these rules of division 

and indexing are required. It is straight forward to consider the modification of this system 
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whereas cell shapes and sections of cells on the plane are included or excluded, bent, joined, 
stretched, rotated, scaled or translated to meet specific spatial requirements such as surface 
modeling or to create tessellations bounded by specific shapes. 

Another flexibility provides for the introduction of one or more new point or cells at a 
5 specific resolution that do not fall on a regular PYXIS location for that resolution. PYXIS 
provides for a cell to be introduced at any unique location and specific resolution. Its ordering 
precedence superceding its neighbours and its behaviors are considered as a parent cell for 
which its centroid location is also a centroid location for lower resolution cells. If two or more 
cells are introduced at any unique location and specific resolution and the boundary of two or 

10 three of the new cells share vertices, such vertices define the location of new child cells and the 
child cells shall be uniquely indexed with reference to its three shared parents, and the behavior 
of these child cells are considered as a parent cell for which their centroid location is not a 
centroid location for any lower resolution cells. An example of this and how it is treated is 
detailed below with the indexing of the first two resolutions of the icosahedron in a discrete 

1 5 global grid system. 

Hierarchal Data Structure 

FIG. 6a and 6b illustrate two additional views of the PYXIS indexing. These figures 
identify a balanced tree hierarchal data structure. To improve the efficiency and speed of 
access for the most probable queries the invention takes advantage of the hierarchical structure 
20 of the dimensions in spatial data. This hierarchical structure reflects the feet that a typical spatial 
query asks about members from the same parent in the hierarchy. Thus, the mapping of the 
present invention provides that records of data with dimension members belonging to the same 
parent are close to each. 

Such data hierarchies in various forms of B-trees and similar data structures are familiar 
25 to those versatile in the art of data management. Spatial data hierarchies are the basis for search 
tree algorithms that allow efficient retrieval of indexed spatial data. 

The identification of nearest neighbour cells is key to such spatial searches. FIG. 7 
illustrates a procedure to identify nearest neighbour cells on the PYXIS Plane (43). The 
following tables are used to provide the foundation for tesseral arithmetic used in Fig. 7. 

30 



14 



WO 2005/013200 



PCT/CA2004/001507 



Table 1 - PYXlndex Class I Matrix 
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Table 2 - PYXlndex Class II Matrix 
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A GBT Gray ordering is exemplary as shown. Addition on Class I (44) and Class II (45) 
5 cell rotation represent respective even and odd indexing digits. The procedure requires use of 
the addition Table (46) whereas the Column Header (47) and Row Header (48) act as the lookup 
values. The Example (49) illustrates the addition of 6 unit vectors to the index value 0600506 
resulting in the nearest neighbour list +1=0600500, +3=0600502, +2=0600050, +6=0600401, 
+4=0605020 and +5=0600504. 

10 Boolean Operations 

Boolean logic plays a central role in spatial data structure algorithms. Boolean algebra is 
useful for performing operations on the attributes (which may be positional or descriptive) 
attached to geographic entities in a Geographic Information System. Boolean Logic is 
especially useful in computing (or modeling) new attributes in topological overlay processing 

15 for both vector and raster based systems, as they can be applied to all data types, be they 
Boolean, Ratio, Interval, Ordinal, or Nominal. A procedure for identifying a PYXIS point as 
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spatially contained, overlapping or excluded, from a specific cell, the spatial data forms of 
Boolean algebra using the logical operators AND, OR, and NOT are outlined below and 
references FIG. 8. 

In the PYXIS Plane (50), any indexed point can be defined as External To (51), 
5 Overlapping (52), or Contained By (53), a Cell (54). Further, the following conditions and 
procedures apply in their determination: 

Condition #1 - A point indexed above a lower resolution containment Cell (54) is 
considered "spatially contained at all resolutions" if its index is nested as a child of the 
containment Cell (54). 
1 0 Test procedure for containment Condition #i 

Given Cell (54) index of resolution n and point index of resolution r, if n < r and digits 1 
to n of Cell (54) is equal to digits 1 to n of the point index then Condition #1 is True. 

Condition #2 - A point indexed above a lower resolution Cell (54) is considered "within 
the spatially overlapping area" if point is found to be "spatially contained at all resolutions" in 
15 one of the six neighbour Cells (56) adjacent to the containment Cell's (54) centroid child Cell 
(55). 

Test procedure for within the spatial overlapping area Condition #2 

Given Cell (54) index of resolution n and point index of resolution r, if n r and if 
digits 1 to n+1 of a neighbour Cell (56) to the centroid child Cell (55) (found using 
20 Addition matrix with unit vectors 1 to 6 from Cell (55)), equals digits 1 to n+1 of the 

point index then Condition #2 is True. 

Condition #3 - A point indexed above a lower resolution containment Cell (54) is 
considered "spatially external at all resolutions" if it is not "within the spatially overlapping 
area" of the containment Cell (54). 
25 Test procedure for within the spatial external area Condition #2 

Given Cell (54) index of resolution n and point of resolution r and Condition #1 is False 
and Condition #2 is False then Condition #3 is True. 

Note that, as a PYXIS indexed point on the PYXIS plane (50) explicitly represents an 
area of influence, a point within the spatially overlapping area (52) may test True under 
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Condition #2 however, if s further determination to be spatially contained, overlapping or 
external to Cell (54) is dependant on the resolution of the point The containment Condition #2 
is therefore resolution dependant Condition #2 only determines that the point is within the area 
of overlapping conditions. To determine if, at the resolution of the point it is contained, 
5 overlapping or external to the Cell (54), further test procedures are required that determine 
where within the overlapping Cell (56) the point lies. 

Within Condition #2, there are two test procedures for cell overlap: Condition #2.1 
whereas the overlap Cell (57) is on a vertex and Condition #23, whereas an overlap Cell (58) is 
on a straight boundary edge. The resolution dependant conditions test procedures are detailed 
10 below and refer to FIG. 8. 

Test procedure for point within a cell overlapping a vertex Condition #2.1 or a cell edge, 
Condition #2.2 

Given point index of resolution r lies somewhere in overlapping n+1 Cell (56), found at 
direction d from Cell (55), test for area that contains point 

15 Fori = 3tor: 

Condition 2.1: 

2.1.1. Point contained in overlapping Area (59) on vertex 

• Check if point contained in n+i centroid child of Cell (57). 

• If point containment True, set i = i + 1 and loop condition #2.1.1 until n+i Cell 
20 index equals point index and therefore "spatially overlapping at this resolution". 

• If False, goto 2.1.2; 

2.1.2. Overlapping Area (60) on edge 

• Use addition matrix to find cell area (60) defined by resolution n+i overlapping at 
edge cells (n+i centroid cells of n+i-1 centroid child of Cell (57) d+2 and d-2 

25 neighbours.). 

• If point containment True, set i = i + 1 and goto condition #2.2 below. 

• If False, goto 2.1.3; 

2.1.3. Overlap Area (61) where all cells contained within the parent Cell (54) 



17 



WO 2005/013200 



PCT/CA2004/001507 



• Use addition matrix to find the 3 resolution n+i contained cells (n+i centroid cell 
of n+i-1 centroid child of Cell (57) d+3 neighbour AND n+i centroid cells of n+i 
centroid child of Cell (57) d+4 and d+5 neighbours.). 

• If point containment True, men point is "spatially contained at all resolutions" of 
5 Cell (54), or; 

• If False, men point is "spatially external at all resolutions" of Cell (54). 

Condition 2.2: 

2.2. 1 . Point contained in overlapping Area (60) on cell edge 

• Check if point contained in n+i centroid child of Cell (58). 

10 .if point containment True, set i = i + 1 and loop condition #2.2.1 until n+i Cell 

index equals point index and therefore "spatially overlapping at this resolution". 

• If False, goto 2.2.2; 

2.2.1. Cells (60) overlapping on edge 

• Use addition matrix to find cell area (60) defined by resolution n+i overlapping at 
15 edge cells (n+i centroid cells of n+i-1 centroid child of Cell (57) d+2 and d-1 

neighbours or d-2 and d+1 if 2.1.1 met in d-2 addition). 

• If point containment True, set i = i + 1 and goto condition #2.2.1. 

• If False, goto 2.2.3; 

2.2.2. Cells (61) contained within the parent Cell (54) 

20 • Use addition matrix to find the 5 resolution n+i contained cells (n+i centroid cells 

of n+i-1 centroid child of Cell (57) d+3 and d+4 neighbour AND n+i centroid 
cells of n+i centroid child of Cell (57) d+3, d+4 and d+5 neighbours.). 

• If point containment True, then point is "spatially contained at all resolutions" of 
Cell (54), or; 

25 • If False, then point is "spatially external at all resolutions" of Cell (54). 
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Spatial Queries 

An embodiment of the invention, as demonstrated above, is a method of selecting the 
spatial data once it resides within the computer device as a hierarchy of indices. Spatial queries 
on the multiresolutional closely packed PYXIS indexing using standard Boolean procedures 
5 familiar to those versatile in the art, provides a means of selecting and sorting spatial features 
and dimensions of their attributes. Examples of practical forms of these selections include: 

• Select all features of a given range of resolutions 

• Select features bound by a given spatial domain (area) 

• Select features bound by a given polygon shape 

10 • Select features from 2 disparate sources A and B, that are in A and in B. (Features 

that exist in both sources) 

• Select features from 2 disparate sources A and B are in A and those in B that are not 
inA. (Add to A all features from B not in A) 

In addition, a multitude of other possible data queries are possible, based on 
15 combinations of resolution, location, cell area, occurrences, attributes, distribution and other 
dimensions of the data. This core functionality provides the basis for further embodiments 
detailed below that include data integration, transmission, visual display, analysis, fusion, and 
modeling. 

Geographical Reference 

20 Another preferred embodiment is a system that utilizes the PYXIS indexing to spatially 

(geospatially) reference (georeference) data located on, around or beneath the surface of earth. 
Methods of projecting points from the faces of a platonic solid to the earth spheroid as a 
geodesic reference are well know in various forms of previous art. A method for establishing 
and referencing multiresolutional closely packed cells on the earth spheroid by projecting the 

25 PYXIS plane from an icosahedron is exemplary of the embodiment and used herein to illustrate 
the capability of the invention, not limiting it to such techniques. 

Prior art described by Sahr, et. al. as the Snyder Equal Area Aperture 3 Hexagon Grid 
(ISEA3H) can be used to orient and generate locations of the closely packed tessellation of the 
icosahedron to the earth spheroid. The PYIXS innovation provides the linear indexing for such 
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a scheme. FIG. 9a shows 5 resolutions of a 3rd aperture hexagon tessellation on a triangular 
face. These cells are indexed in a manner as described above for linear storage and hierarchal 
reference. FIG 9b shows an icosahedron solid featuring 12 Vertices and 20 triangular Faces. 

Reference to the vertices requires a modification of the PYXIS innovation such that the 
5 shape of the original vertex cells and all subsequent cells that are created at this same centroid 
(always at the icosahedrons vertex) will be a pentagonal as illustrated in FIG. 9c. Further, the 
cells are laid onto the five feces related to each vertex of the icosahedron. The indices are 
modified accordingly. 

FIG. 9d shows 6 resolutions of the modified PYXIS plane oriented on an unfolded (2D) 

10 icosahedron suitable for simple transformations. FIG. 9e shows this same plane reference 
projected to a sphere using John Snyder's modified Lambert Azimuthal Equal Area Projection. 

In addition to the typical PYXIS tree hierarchy developed above and shown in FIG. 6a 
and 6b, the indexing hierarchy in the georeference to the icosahedron must be modified to 
include two additional levels, one to index the vertices and a second to index the points on the 

15 faces. Exemplary of a simple ordering for these points, but not limited to these labels, would 
have the vertices labeled 1 through 12 and the faces labeled A to J and Q to Z. This face labels 
are further referenced to their three parents to provide georeference and order. As an example 
child cell A is related to parents 1, 3 and 5. Hierarchal precedence is given to the modified cells 
deemed to be resolution 1 created at the vertex of the icosahedron where the vertices of this 

20 modified cell are located at each of the five connected icosahedron faces. Reference at 
subsequent resolutions returns to the regular PYXIS indices. 

In this way, the PYXIS innovation is used to provide a closely packed, uniformly 
adjacent, multiple resolution, overlapping spatial data discretization and ordering that can be 
used to georeference the earth spheroid and allow spatial observations of the earth to be 

25 transposed from planar space to a one-dimensional space efficient for use in digital computer 
memory and storage. 

A Discrete Global Grid System 
System Environment 

FIG. 10 illustrates the components of a Discrete Global Grid System (63) as a further 
30 embodiment of the invention. The PYXIS innovation provides the combined utility of a closely 
packed, multiresolutional cell index with a georeference to the earth spheroid that has 

20 
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hierarchical features suitable for storage and selection by structured query. The PYXIS 
innovation is the basic framework for components of the system: localized and global spatial 
database (64) (storage (70) and query (71)), retrieval (65) (conversion (72), geo-referencing (73) 
and integration (74)), data discovery and file sharing (66) (compression (75) and transmission 
5 (76)), data transformation systems (67) (visual display (77), analysis (78), fusion (79), and 
modeling (80)). Creation of imagery and data can be completed by raster and vector image 
sensors, creators and editors (68) or as a product of the Discrete Global Grid System (63). Peer- 
to-peer networking of PYXIS data discovery components allow data to be advertised and shared 
between remote computers (69). 
10 There are many specific methods known within the art of spatial data storage, retrieval, 

transmission, visual display, integration, fusion, analysis, and modeling suitable for use with the 
PYXIS innovation. Further, there are other functions that build on the PYXIS innovation not 
specified. The following is summarily described functionality and it should not be interpreted 
as encompassing all functional applications of the PYXIS system. 

1 5 Spatial Data Retrieval 

The Spatial Data Retrieval System allows images and data of various formats, sources, 
types and georeference to be converted to PYXIS geospatial indices, through binning and 
gridding algorithms, suitable for storage. The critical component to this procedure is the 
capability to convert a Point (81) with a rectangular coordinate to a cell with a PYXIS index. 
20 This can be accomplished by determining the index of the PYXIS cell that contains the point. 
The point's area of influence provides a basis for assigning a cell resolution. 

With reference to FIG. 11, the following steps provide an exemplary method of 
determining if a rectangular point is contained in a specific PYXIS cell: 

1. Align the rectangular Plane (82) with the PYXIS Plane (83) such the rectangular 
25 coordinate of at least one cell centroid is known. Let <|> be the angle between the PYXIS 

Plane (83) and the rectangular Plane (82) and a is 30° for cell resolutions with vertex 
oriented at the top and 60° for cell resolutions with cell edge oriented at the top of the 
PYXIS plane. 

2. Determine the intercepts of a Y axis with two lines from Point (81), one at an angle -(a 
30 + <j>), referred to as -YintpomT and the other at a - cj>, referred to as YintpoiNT. 
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3. Determine the intercepts of the same Y axis with two lines from the centroid of Cell 
(84), one at an angle -(a + <|>), referred to as -YintcELL and the other at a - <j>, referred to 
as YintcEix- 

4. Where: 

5 a is die length of one cell side 

b is -YintpoiNT minus -YintcELL 
c is YintpoiNT minus YintcELL 
k = a/cos(a + <j>) 
j=a/ cos(a - <J>) 
10 Then Point (81) is contained in Cell (84) if: 

-j > b > j AND -k > c> k AND b - c> a. 
Note that in the example shown, Point (81) is found not contained in Cell (84). 
Data Visualization 

Data is indexed and can be stored in a database management system or PYXIS file. As 
15 described above, spatial features can be efficiently selected by their spatial relationships, 
features and dimensions to complete further transformations of the data. 

Digital Globe Visualization is an embodiment achieved by transforming selected 
features from a stored PYXIS index to a rendered object or graphic. Selected feature geometry 
and attribute values are sent to a rendering pipeline. The application returns to a computer 

20 visualization device a representation of the spatially organized data associated with a spatial 
area and range of resolutions in the form of a whole or partial rendered image of the geodesic 
globe. The resultant Digital Globe is capable of being rotated, translated, and scaled to produce 
zoom, pan, yaw, and rotation functions. Exemplary data types eligible for display include 
PYXIS coded raster images, vector features, map annotations, texture bitmaps and terrain 

25 surface and 3 dimensional models. 

Geospatial Data Discovery 

Another preferred embodiment of the PYXIS Discrete Global Grid System is an Online 
Spatial Data Discovery and File Sharing Server. With the addition of georeferenced PYXIS 
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index to multidimensional attributed, fielded, and textural or image data, information can be 
discovered based on its content and geospatial position. Prior art forms using standards and 
protocols for data exchange and display, like HTML, XML, SOAP, and other web services or 
peer-to-peer file sharing provide the basis for spatial data publishing and discovery. 

5 Further, the file of data can then be transmitted from the server on demand in its 

entirety, progressively or continuously. PYXIS code data resides locally and/or remotely and 
can appear on demand, progressively transmitted or updated continuously over a network 
environment 

As an example, a schoolteacher, who wants to show her class of students a computer 
10 view of the globe, searches, discovers, and then accesses a remote site that serves a vector file of 
world political boundaries; gained at another server site, she overlays the globe with a recent 
mosaic of satellite images. Setting resolutions for display, she downloads the image on- 
demand. The large satellite image is served ''progressively", based on the ideal visual resolution 
and the extent of the view on die computer monitor. As the teacher zooms into a specific area, 
15 say downtown Ottawa, the extent of the satellite image is decreased and the resolution increases. 
At some point a high-resolution digital ortho photo of Ottawa replaces the lower resolution 
satellite image. Level MapO vector file, a 1:100,000 digital topographic map of the world, 
replaces the boundary file. This type of transmission is suitable as streaming data. 

Another form of file transfer, continuous transmission, is helpful when sensor data is 
20 used to update a real-time event While observing a road intersection in downtown Ottawa, a 
round symbol is changed on the globe from red to green when a traffic signal changes, or a 
temperature sensor provides for changes in annotation from 22C to 21C. 

Another example of a transformation embodiment is the basis for file transmission. A 
feature of the PYXIS innovation is the compactness of the index. First level compression of a 
25 series of indices for transmission is accomplished by noting that cells of different resolutions but 
in the same area retain the same parent. By ordering the cells in their hierarchy from lowest 
resolution to highest the parent cell of any index need only be sent once if sent in order. As an 
example and not limited to this technique, given the following indices: 

010203, 010204, 01020401, 01020402, 01020405, 01020500, 01020501, 01020502, 
30 01020503, 010306, 0105 
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The following digits can be transmitted without losing data where 0 indicates branching back up 
to the last parent Oast 0): 

010203, 4, 01, 2, 5, 0, 500, 1, 2, 3, 0, 0, 036, 0, 0, 05 

Further compression encoding, whether lossy or lossless, can be performed on PYXIS 
5 indices by methods known by those versatile in the art. Further, pre-queried items can serve as 
a dimension in a spatial OPAP data cube or the data set for a particular CODEC. 

Files encoded in this way can be transmitted in their entirety from a server to a client A 
PYXIS index attached to a device location, like an IP address, or embedded in a HTML file, or 
as an XML type tag, provides a geo-reference to the data contained in the file. In this way, the 
10 digital representation of such data can be referenced to me earth, providing a spatial abstraction, 
the globe, as a Spatial Data Browser (or Viewer or Media Player) is an embodiment of the 
invention. Data stored in PYXIS format are identified in an on-line environment as a search 
query, displaying at an automated or manually set resolution, a pictographic symbol at the cell 
location on the Digital Globe. Users can select this symbol with a cursor, activating further 
1 5 software instructions or viewing detailed data. 

Or, the parts of me index and attributes may be transferred from a server to a client in 
response to a query that bounds the area and range or resolutions described above as a query of 
the PYXIS data. This "progressive transmission" would be activated manually or could be 
automated based on properties of the clients viewing area. As a user "zooms" from a larger to 
20 smaller area, the system returns data further down the resolution of the spatial hierarchy, less 
area, more detailed. 

Or, the data could be continuously transmitted, wherein the client is continuously 
checking a server for changes to the attributes of a specific cell or group of cells. An example 
may be an updated RGB colour value from a georeferenced PYXIS indexed video, photometric 
25 or satellite image or any other sensor device that includes PYXIS cell reference. 

Or, the data can be queried and built as streaming data providing suitable options for 
real-time applications. 
Spatial Data Analysis 

An embodiment of the system is Spatial Data Analyzer. Data Analysis Transformations 
30 in PYXIS are achieved by Boolean query, performing operations on the attributes (which may 
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be positional or descriptive) to select, build and enhance geography. Spatial analysis with the 
PYXIS DGGS relies heavily on the geo-synchronization of the gridded data. As an example, 
using the temporal dimension of census data, a comparison of population of a city in the year 
2002 with its population in 1950 is strictly correct if the boundary of the city remained 
5 unchanged. Gridded data analysis provides an efficient method of modifying that query to a 
comparison of the population of a specific set of cells in 2002 to those same cells in 1950. It is 
thus that PYXIS can perform analysis. 

Data Fusion Transformations in PYXIS provide real interpolation of data. The fusion 
functionality allows the user to select features for 2 or more sources of data, similar to spatial 
10 analysis, and then create new attributes through statistical methods known to those versatile in 
the art of data fusion. Data Integration Transformations, which are similar to data fusion, are 
achieved by selecting data that are included in a polygon and performing geometric transforms 
(move, scale, rotate) and superimposition (move data above or below, cut and paste, overlay and 
mosaicking). 

15 Data Modeling Transformations in PYXIS provides real extrapolation of data. Many 

studies of dynamic environmental systems can be semi-deterministically modeled with access to 
current accurate data referenced to a common grid along with formulation of cellular behavior. 
Modeling methods such as finite difference and cellular automata rely on a defined mesh and or 
tessellation to frame the calculations. 
20 An embodiment of the invention is a Geospatial Model Building system, which allows 

the closely packed, uniformly adjacent, multiresolutional, overlapping spatial data ordering to 
be used as a mesh and grid for the construction of stochastic and deterministic simulation of 
dynamic earth events. Users can access on-demand in a peer-to-peer environment a multitude 
of temporal geospatial data at each cell and extract and utilize this spatial data. Users can access 

25 properties of the grid and the data that references it. Formulas, with constants and variables 
with spatial and temporal values assigned from PYXIS data. Rules of behavior can be defined 
which allow extrapolation and predictive modeling, simulating future cell conditions. 

As an example, in a PYXIS system, data pertinent to rainfall runoff modeling is 
discovered, extracted and transmitted from a variety of sources to the model builder including 

30 current rainfall conditions, temperature, vegetative cover, soil types, terrain, etc. The user 
builds formulas using these data fields and other coefficients as variables for calculating water 
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budget, stage, storage and routing conditions. The simulation of a rainfall event results in the 
increase of simulated water flow over time in a downstream cell, emulating a rainfall/runoff 
hydrograph. 

While the present invention has been described and illustrates herein with reference to 
5 the preferred embodiment thereof it will be understood by those skilled in the art that various 
changes in form and details maybe made therein without departing from the spirit and scope of 
the invention. 

It is to be understood that the embodiments and variations shown and described herein 
are merely illustrative of the principles of this invention and the various modifications may be 
10 implemented by those skilled in the art without departure from the scope and spirit of the 
invention 
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